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ABSTRACT 

We have investigated the evaporation of close-in exoplanets irradiated by ionizing photons. We 
find that the properties of the flow are controlled by the ratio of the recombination time to the flow 
time-scale. When the recombination time-scale is short compared to the flow time-scale the the flow is 
in approximate local ionization equilibrium with a thin ionization front, where the photon mean free 
path is short compared to flow scale. In this “recombination limited” flow the mass-loss scales roughly 
with the square root of the incident flux. When the recombination time is long compared to the flow 
time-scale the ionization front becomes thick and encompasses the entire flow, with the mass-loss 
rate scaling linearly with flux. If the planet’s potential is deep the flow is approximately “energy- 
limited”; however, if the planet’s potential is shallow we identify a new limiting mass-loss regime, 
which we term “photon-limited”. In this scenario the mass-loss rate is purely limited by the incoming 
flux of ionizing photons. We have developed a new numerical approach that takes into account the 
frequency dependence of the incoming ionizing spectrum and performed a large suite of ID simulations 
to characterise UV driven mass-loss around low mass planets. We find the flow is “recombination- 
limited” at high fluxes but becomes “energy-limited” at low fluxes; however, the transition is broad 
occurring over several order of magnitude in flux. Finally, we point out the transitions between the 
different flow types does not occur at a single flux value, but depends on the planet’s properties, with 
higher mass planets becoming “energy-limited” at lower fluxes. 

Subject headings: hydrodynamics, radiative transfer, planets and satellites: general, ultraviolet: stars 


1. INTRODUCTION 

Observational campaigns to detect exoplanets have 
been incredibly successful, having found thousands of 
possible exoplanets (e.g. Batalha et al. 2013). How¬ 
ever, what’s surprising is that many of the exoplanets are 
found to be incredibly close to their parent stars with sep¬ 
arations < 0.1 AU being common. Furthermore, while 
the original detected planets at small separations were 
massive (of order a Jupiter mass and higher) - so called 
“hot jupiters” - these are relatively rare (Fressin et al. 
2013). With the advent of the Kepler mission the major¬ 
ity of planets found close to their parent star are small 
radius, low mass objects, with the majority of stars con¬ 
taining at least one small (< 4 R^) planet with a period 
shorter than Mercury’s (Fressin et al. 2013). This in¬ 
dicates that the formation of close-in small planets is a 
dominant mode of planet formation. 

At such close separations, the environment is very dif¬ 
ferent to planets in our own solar system and irradiation 
from the central star can play a very important role, par¬ 
ticularly when the star is young (e.g. Owen & Wu, 2013) 
or during the late stages of stellar evolution (e.g. Bear & 
Soker, 2011; 2015). The star’s bolometric luminosity can 
heat-up the planet’s atmosphere leading to an increased 
scale height, as well as preventing the planet from cool¬ 
ing (Guillot et al. 1996, Burrows et al. 2000, Baraffe 
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et al. 2003). Furthermore, the high energy ionizing lu¬ 
minosity from the star can heat-up the upper layers to 
temperatures of order the escape temperature, leading 
to mass-loss (e.g. Lammer et al. 2003, Tian et al. 2005, 
Owen & Jackson, 2012). 

When the star is relatively young (< several 100 Myr 
old) the ionizing luminosity can represent a significant 
fraction of the total output of the star ~ 10~ 3 L* (Jack- 
son et al. 2012). Integrating the total received ioniz¬ 
ing luminosity over Gyr time-scales and comparing it 
to the binding energy of close in planets it was realised 
that planetary evaporation could play a significant evo¬ 
lutionary role and actually remove a large fraction of the 
planet’s atmosphere (Lammer et al. 2003, Jackson et al. 
2012), especially at low planet masses (Owen & Jackson 
2012, Lopez et al. 2012, Owen & Wu 2013, Lopez & 
Fortney 2013). 

Due to the importance of the evolutionary effects of 
planetary evaporation one requires accurate mass-loss 
rates. Commonly the approach is to turn a certain frac¬ 
tion of the received ionizing luminosity into work/heating 
and associated mass-loss, either globally (Lammer et al. 
2003, Erkaev et al. 2007) and obtain the well known 
“energy-limited” mass-loss rate; or locally in simplified 
simulations (e.g. Yelle et al. 2004, Tian et al. 2005, Lam¬ 
mer et al. 2013). Such an approach implicitly assumes 
that PdV work dominates the energy loss and radiative 
cooling or ionization physics are unimportant. 

More detailed studies have since been performed that 
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take into consideration some of the details of radiative 
transfer and ionization chemistry (Koskinen et al. 2007, 
2010, 2014; Chadney et al. 2014); however, these mod¬ 
els typically concentrate on a handful of observed hot 
jupiters. Owen & Jackson (2012) considered the role X- 
rays play and showed heating from secondary ionization 
and line cooling were important when the sonic point 
occurred in an area of the flow dominated by X-ray 
heating and argued that in this case the flow was never 
“energy-limited”. Owen & Jackson (2012) also presented 
a schematic picture of when UV or X-ray heating dom¬ 
inates and argued that X-rays would dominant for the 
highest fluxes around low-mass planets. 

Furthermore, Murray-Clay et al. (2009) considered the 
evaporation for a hot jupiter with mass 0.7 Mj and radius 
1.4 Rj with a pure UV radiation field at various fluxes. 
At high fluxes > 10 4 erg s _1 cm -2 , they found the flow to 
be in approximate radiation-recombination equilibrium, 
in this case the flow behaved like a HII region, was fully 
ionized and approximately isothermal with a tempera¬ 
ture of 10 4 K and the cooling was dominated by Lyman- 
a radiation. As the density in any HII region scales with 
as </q/ 2 (where Jo is the ionizing flux in photons per unit 
time), then the mass-loss rate in this regime scale ap¬ 
proximately as J 0 1/2 . At fluxes < 10 4 erg s _1 cm -2 the 
density in the flow was so low that recombination were 
too slow compared to the flow time and the mass loss 
rate then scaled linearly with flux. In this case ioniza¬ 
tions were balanced approximately by advection, which 
Murray-Clay et al. 2009 identified as an “energy-limited” 
flow. 

The physical understanding of the transition for 
“recombination-limited” to an “energy-limited” flow be¬ 
ing that once the density drops in the flow such that the 
recombination time is long compared to the flow time 
then it can no-longer be in radiative-recombination equi¬ 
librium (e.g. Soker, 1999). This obviously depends on 
planet properties, through the flow time (e.g. planet 
radius and mass); and as such the transition flux be¬ 
low which the flow cannot be in radiative-recombination 
equilibrium is not going to be a unique value. However, 
several studies of exoplanet evolution (e.g. Jin et al. 
2014) take the results of Murray-Clay et al. (2009) for 
their single planet study and apply it to all planets which 
is likely to lead to widely inaccurate adopted mass-loss 
rates. 

In this paper we aim to study the properties of UV 
driven evaporation across the range of observed exoplan¬ 
ets, particularity at low masses which has as yet been un¬ 
explored. We aim to put forward a simple physical under¬ 
standing of the various evaporation regimes and identify 
three: “recombination-limited”, “energy-limited” and a 
new regime “photon-limited”. We present an analytic 
picture in Section 2, perform a numerical study with 
a new radiation-hydrodynamic scheme in Section 3. In 
Section 4 we discuss our results and summarise in Sec¬ 
tion 5. 

2. ANALYTIC OVERVIEW 

The evaporative flow problem is complicated by the 
fact that in general the flow may not be in local thermo¬ 
dynamic or ionization equilibrium, particularly at low 
fluxes (e.g. Murray-Clay et al. 2009). However, in the 


case of irradiation by a spectrum dominated by Hydro¬ 
gen ionizing photons (EUV) we can make progress by 
considering ionization balance (e.g. Bertoldi & Mckee, 
1990). We consider a planet with mass M p and radius 
R p close to its central star at a separation a. Noting 
that R p /a <C 1, then the EUV flux is constant on all 
scales of interest. We consider a flow that is dived into 
three distinct regions: (i) bolometrically heated planet’s 
atmosphere; where the temperature is set by the flux of 
optical photons and the optical depth to EUV photons 
is large such that the ionization fraction (A) is very low. 

(ii) an “ionization front” where the photo-ionization rate 
greatly exceeds the recombination rate and the ionization 
fraction increases as one moves away from the planet. 

(iii) a fully ionized region where photo-ionization is bal¬ 
anced by recombinations and the ionization fraction is 
unit50 

Since the scale height of the planet’s bolometrically 
heated atmosphere must be <C R for a bound planet 
then the transition from the bolometrically heated at¬ 
mosphere (region i) to the ionization front (region ii) is 
by construction sharp. In the picture we develop below 
we take it to occur at an optical depth teuv = 1 to the 
EUV photons and call this radius the planet’s radius. 
However, unlike many problems in astrophysics involv¬ 
ing EUV radiation, as we shall see we cannot in general 
assume that the ionization front is thin compared to the 
flow scale (or even the planet’s scale) and in many cases 
we find a transition from region ii to region iii occurs at 
a very large radius. In particular, if it occurs at a radius 
larger than the sonic radius then it is irrelevant in de¬ 
termining the mass-loss and we assume the flow to only 
contain regions i and ii. 

We choose to work in one-dimensional spherically sym¬ 
metric picture along a streamline connecting the star and 
planet (see Stone & Proga, 2009 and Owen & Adams, 
2014 for the role of multi-dimensional flows). For the 
sake of clarity we neglect the contribution from the tidal 
gravity due to the star. For a steady-state, spherically 
symmetric flow the governing equations of hydrodynam¬ 
ics are: 
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= q 
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where p is the 

mean-molecular weight given by 

M = 


mn/{ 1 + A), with A' the ionization fraction; G is the 
gravitational constant; n is the number density of Hydro¬ 
gen nuclei; p = np the mass density; u the gas velocity; 
P = ki,nT the gas pressure, where T is the tempera¬ 
ture; E is the total specific energy density and q is the 
net volumetric heating rate, which is detailed in Section 
3.1. Furthermore steady state ionization balance requires 
(e.g. Bertoldi & McKee 1990, their Equation 2.1): 

rRp 

Jo=J r w- / dRa R X 2 n 2 (4) 

J OO 

2 Here we consider the case of a pure Hydrogen atmosphere, with 
the inclusion of other elements the ionization fraction is generally 
> 1 
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where an is the recombination rate and Jp is the flux of 
ionizing photons into the ionization front. Since recombi¬ 
nation are negligible inside the ionization front then Jp 
must balance the flux of neutral particles into the ion¬ 
ization front (Churchwcll et al. 1987, Bertoldi & McKee 
1990). Namely: 

Jif = n(Rp)u(R p ) (5) 


Note since the radiation held is plane parallel, Equa¬ 
tion [4] is written in plane parallel form. Since the flow 
is approximately spherical we must be careful combin¬ 
ing Equation [4] with the hydrodynamic terms, something 
we discuss further in Section 2.1.1. Specifically, we can¬ 
not trivially integrate Equation [5] over the surface of the 
planet by multiplying by AttR^ to find the mass-loss rate, 
since the total number of ionizing photons entering the 
ionizing front does not necessarily equal the total num¬ 
ber reaching the base of the atmosphere, as some may 
pass through the how without interacting with the it. 
Equation [3] can be recast into a limiting mass-loss, the 
“energy-limited” mass-loss rate where uncertainties in 
heating and cooling are encapsulated into an efficiency 
parameter 77 which can only be determined through nu¬ 
merical calculations with values typically in the range 
0.01-0.3 (e.g. Murray-Clay et al. 2009; Owen & Jackson, 
2012 ). 


M < i) 


7 xR? p J 0 hv h 
AGMp 


( 6 ) 


where hVh is the mean photon energy that goes into heat¬ 
ing (rather than ionization). The origin of the factor 4 
in the denominator comes from averaging the received 
flux by a nR^ disc over the dnRp surface area of the 
planet. We choose to this form as it is the commonly 
used form in the exoplanet literature and allows compar¬ 
isons of the efficiency ( 77 ) values, although we note some 
authors do not perform such an averaging (e.g. Murray- 
Clay et al. 2009). For a flow that can be considered close 
to being energy-limited requires that q is predominately 
controlled by heating and radiative cooling is a small 
contribution. Therefore, energy losses are dominated by 
PdV work or advection. The main contribution to the 
radiative cooling rate A is Lya radiation with cooling 
rate (Black 1981): 


A = 7.5 x 10 19 erg cm 3 s 1 


n e no exp 


-118348 K 
T 


(7) 

The ionization state of the gas is controlled by the 
recombination rate; adopting the on-the-spot approx¬ 
imation such that aR = as where olb is the case- 
B recombination coefficient given by as = 2.6 x 
10 -13 cm 3 s _ 1 (T/10 4 K ) -0 - 7 so that the recombination 
rate Ai on is: 

Aion = &B n e (8) 


where n e is the electron density. Noting the ratio of the 
cooling time-scale ( t c ) to the recombination time-scale 
(t r ) is given by: 


tc 


LbT Ai on 

A 


(9) 


and is purely a function of temperature and ionization 
fraction. When the gas is close to ~ 10 4 K and roughly in 


ionization-recombination equilibrium then t c /t r ss 0.05. 
However, if the flow were no longer in local ionization 
equilibrium and the temperature falls below 10 4 K. The 
temperature falls because the heating rate and ionization 
rate scale identically with the ionizing flux; therefore, 
since the ionizing photons are unable to obtain ioniza¬ 
tion equilibrium they are unable to heat the gas to 10 4 K 
temperature (i.e. the ratio of the heating to ionization 
time-scale is a constant). Thus, PdV cooling becomes 
the dominant cooling mechanism (e.g. Murray-Clay et 
al. 2009), so we find t c /t r ~ exp(—10 5 K /T). Thus, 
t c quickly becomes longer than the recombination time 
when the flow is no longer in local ionization equilib¬ 
rium. Since local ionization balance breaks down when 
when the recombination time becomes longer than the 
flow time-scale ffl ow , then one would expect that in the 
majority of the cases t c > ffl ow when T < 10 4 K and the 
flow will be close to “energy-limited” (c.f. Soker 1999, 
Bear & Soker 2011). This does not mean that Equations 
will provide an accurate description of the mass-loss rare 
over the range of parameters of interest as 77 may still 
vary considerably due to hydrodynamic effects. 

2 . 1 . Limiting cases 

Inspection of Equation |i]&[ 6 ] indicate that three phys¬ 
ical processes may limit the mass-loss. Either the num¬ 
ber of incoming photons, recombinations or energy con¬ 
siderations. Equating Jq with the advection term in 
Equation [4] indicates that the flow will be limited by 
the number of incoming ionization photons — a pro¬ 
cesses we term “photon-limited”. Equating Jq with 
the recombination term in Equation Id] results in a den¬ 
sity in the ionization front that is determined only by 
radiation-recombination balance — a processes we will 
term “recombination-limited”. Finally, insisting that the 
flow satisfies Equation [ 6 ] will result in a flow that is lim¬ 
ited by the available energy, this has become known com¬ 
monly as “energy-limited” evaporation. In this case the 
mass-loss rate is described by taking Equation [ 6 ] as an 
equality where 77 is yet to be determined. 

2 . 1 . 1 . Photon limited flows 

This corresponds to the case where the ionization front 
encompasses the entire flow (including the sonic point) 
and as such recombinations are completely negligible in 
the flow so that Jif = Jq = n(R p )u{R p ). Additionally 
the pull of gravity is weak and the inequality in Equa¬ 
tion 6 is easily satisfied. Note we are implicitly assuming 
that the gas temperature is high enough that the gas is 
able to escape freely. . Therefore, we are limited by the 
number of incoming ionizing photons: namely we can 
only loose one atom per incoming photon. Note we are 
implicitly assuming that the gas temperature reached in 
the ionization front is high enough that escape from the 
planet’s potential is possible; if not the flow becomes 
“energy-limited” as we will discuss in Section 2. How¬ 
ever, here we are complicated by geometric effects. For 
uniform spherical irradiation then the total number of 
absorbed photons is the area out to the ionization front, 
not the planet’s radius itself. However, planet evapora¬ 
tion does not occur with spherical irradiation, but rather 
plane parallel. Furthermore, since the ionization front 
is marginally optically thin ionizing photons, then any 
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photons with an impact parameter much larger than the 
planets radius, but smaller than the ionization front ra¬ 
dius are likely to pass straight through the flow. There¬ 
fore, we can approximate the effective absorbing area of 
the planet as ~ nR p . Therefore, in order to obtain the 
photon limited mass-loss rate we drop the recombination 
integral in Equation [4] and equating J$ and the advection 
allows us to calculate the mass-loss rate approximately 
as: 

M PL « tt R 2 p m H J 0 (10) 

We note that we do not need to make assumptions about 
the energy considerations here as we do not need to know 
the gas temperature to find the mass-loss rates. Such a 
mass-loss limit has not previously been discussed in the 
context of planet evaporation. Thus, it will corresponds 
to cases where planetary gravity is not strong and the 
flux is low, so namely low-mass planets at large separa¬ 
tions and early times. Unfortunately, due to the fact any 
ID model essentially assumes spherically symmetric irra¬ 
diation any photon limited mass-loss rate is likely to be 
much higher than this value (since the effective absorb¬ 
ing area in such a calculation is 7ri? 2 F and Rif > R P 
often with R\p R p ). Therefore, in our numerical 
calculations presented in Section 3 we will not be able 
to study this limit of planet evaporation. Any numeri¬ 
cal constraints on “photon-limited” evaporation must be 
obtained with 2/3D numerical calculations. Although, 
Equation [lO] should roughly represent an appropriate es¬ 
timate for the mass-loss rate, and is similar to the 1+1D 
numerical calculations of Bertoldi & McKee (1990). 

2.1.2. Recombination limited flows 


Writing the term in parenthesis in Equations El and[l2] 
as a scale height H then we crudely set: 


H = min 


R p 

( c 2 R 2 p V 

3 ’ 

\2 GM P J 


(13) 


Therefore, we can write the base density as: 


ni = v5? (14) 

Now we can write the recombination-limited mass-loss 
rate, noting in fully ionized flow p, = m R /2 as: 

Mrl = 2nR 2 p m H n b c s M b (15) 

where A4 b is the Mach velocity at the base of the 
flow. For an isothermal flow M b is given by the well 
known Parker wind solution (Parker, 1958, Cranmer 
2004) where: 


M b = J-W 0 


Rr, 


exp 3 —— 1 


R. 


4 R„ 


Rr, 


(16) 


with W(x) the Lambert W function and R s = GM p /2c 2 , 
if Rp/R s < 1 and M b = 1 otherwise. For tightly bound 
planets Equation 16 has an exponential fall off of the 
form: 



(17) 


Therefore for R p /R s close to unity M RL oc i? 3 / 2 ; however 
for Rp/R s -C 1 it has a steep (exponential) dependence 
on mass and radius. 


“Recombination-limited” flows are perhaps the easiest 
to understand as in this case the radiative transfer and 
hydrodynamics problems are essentially de-coupled and 
the basics in the context of planetary evaporation were 
discussed in Murray-Clay et al. (2009). When the recom¬ 
bination time is short compared to the flow time-scale the 
ionization front becomes thin compared to the flow scale 
and essentially the transition from region i to region iii 
can be considered infinitesimally thin. Furthermore, the 
gas is tliermostated to approximately 10 4 K and approx¬ 
imately isothermal. For cases when the planet’s gravity 
is strong enough that the scale height in 10 4 K is small 
compared to the planet’s radius then one can approxi¬ 
mate the density structure close to the planet as hydro¬ 
static and the integral in Equation [4] becomes (see also 
Murray-Clay et al. 2009): 

( 11 ) 

where n b is the density at the ionization front. However, 
in the case that the planet’s gravity is weak then the 
density profile will be roughly that of a constant velocity 
outflow and will fall of as r~ 2 . Therefore, in this case 
the integral in Equation E] becomes (see Johnstone et al. 
1998): 

J dRa R X 2 n 2 « a B n b ( 12 ) 


2 .2. Picture of Evaporation across parameter space 

For fixed mass, since Mel oc r]R p 3 and M PL oc R 2 , 
(i.e. they intercept once) then we note that there are 
two possible cases. Either the flow is “recombination- 
limited” at all possible radii (high masses and high 
fluxes) or the flow is “recombination-limited” at both 

large radii (where M R l oc R p ^ 2 ) and at small radii 
where it has an steep dependence on radius. Further¬ 
more, since “photon-limited” mass-loss has a shallower 
dependence on radius than “energy-limited” then photon 
limited flows will dominate at larger planetary radii and 
lower masses than “energy-limited”. The transition be¬ 
tween “photon-limited” and “energy-limited” flows can 
be found by taking Equa tion [6| as an equality and equat¬ 
ing it with Equation |10| to find: 

Mp = V (dGm J Rp (18) 

Note this expression indicates that when the photon en¬ 
ergy is much larger that the depth of the gravitational 
potential, then the photon limit is the relevant one at low 
fluxes (i.e. we can assume heating is sufficient to allow 
the atmosphere to escape as we adopt in Section 2.1.1). 
Whereas, when the photon energy is much less than the 
depth of gravitational potential, the energy limit is the 
relevant one at low fluxes as one cannot assume the gas 
can escape freely. 
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One can find the transition between “recombination- 
limited” and both “photon-limited” and “energy- 
limited” by equating Equation 15 with Equations 10 & 
[6] respectively to find: 

-4 
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exp 
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JoubH 
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(19) 


for the transition between “recombination-limited” and 
“photon-limited” and 
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for the the transition between “recombination-limited” 
and “energy-limited” where we note the last term in 
brackets is the rat io o f the photon energy to binding 
energy. Equations [19] & [20] posses either two solution 
or no solutions and cannot be presented in closed form. 
Therefore, we proceed graphically. 

In computing the various mass-loss regimes we assume 
that the mass-loss efficiency is constant with planet ra¬ 
dius and mass. As lowering the efficiency reduces the 
energy limited mass-loss rate this makes this limit more 
restrictive and increases the range of parameters over 
which the flow can be characterised by energy-limited 
evaporation. For our fiducial case we take hi>h = 6.4 eV, 
rj = 0.1 and use an EUV flux of 3 x 10 4 erg s” 1 cm -2 and 
the various mass-loss regimes are shown in top left panel 
of Figure [T] In our panels of Figure [T] we also show a 
recent census of the exoplanet population as points. 

This moderate flux corresponds roughly to a separation 
of 0.2 AU around a young star and 0.04 AU for a 1 Gyr 
old star (Jackson et al. 2012). In this example we see 
that jupiter mass and above planets are evaporating in 
the recombination limit, whereas most sub-jupiter mass 
planets are losing mass in an energy limited fashion. For 
very low mass planets (~M$) we find that evaporation is 
photon limited. We note that while most planets evolve 
to the left in the diagram as they shrink this is only a 
minor adjustment for massive planets. However, planets 
with envelope mass-fractions of order 10s of percent can 
be much larger at early times, and can often have radii 
5-10 times large than at Gyr ages when they are observed 
today (Lopez & Fortney 2013; Owen & Wu 2013). 

In the bottom left and top right panels of Figure [T] we 
explore the role incident flux plays, with fluxes of 3 x 10 3 
and 3 x 10 5 erg s _1 cm -2 respectively. As one expects at 
lower fluxes the parameter space dominated by energy- 
limited and photon-limited flows, with the majority of 
planets loosing mass in an energy limited fashion and a 
small minority in a photon limited picture. While mas¬ 
sive (> Mj) planets still remain in the recombination 
limit. 

Alternatively at high fluxes no planets are in the 
photon-limited picture and the energy limited flow occu¬ 
pies a small fraction of parameter space and most planets 
are evaporating in the recombination limited scenario. 
Since at this flux the majority of planets that sit in the 
energy limited-region are low mass. In order to experi¬ 
ence fluxes of this magnitude requires the planets to been 
irradiated at early time (< 100 Myr - Jackson et al. 2012, 


Owen & Wu 2013) and as such the planets plotted here 
would have been considerably larger and most would sit 
to the right of the energy-limited region. 

Finally, we investigate how the mass-loss efficiency pa¬ 
rameter in the energy-limited flows affects the param¬ 
eter space. In the bottom right panel of Figure |T| we 
show our fiducial model, but with an efficiency parame¬ 
ter of ij = 0.3 corresponding to the upper range found by 
Owen & Jackson (2012). As expected increasing the effi¬ 
ciency parameter reduces the importance of the energy- 
limit and the majority of low mass planets are now in 
the photon-limit. Whereas energy-limited evaporation is 
constrained to intermediate mass-planets. 

The picture described here is consistent with the study 
of Murray-Clay et al. (2009) who studied a 0.7 Mj 
planet with a radius of 10 10 cm as a range of fluxes and 
found that the flow transitions from recombination lim¬ 
ited at high fluxes to energy limited at low fluxes around 
~ 10 4 erg s _1 cm -2 , consistent with our calculations 
that predict a transition flux of ~ 2 x 10 4 erg s _1 cm -2 . 

Now we can build up a picture for the mass-loss his¬ 
tories of both low mass planets (< 2 OM 0 ) and high 
mass planets (~Mj). Close-in high mass planets begin 
evaporating in the recombination limited regime, as they 
evolve and the EUV flux falls with time they may tran¬ 
sition to energy limited at some low flux, for typical pa¬ 
rameters this is about 10 4 erg s^ 1 cm -2 ; however, we 
note that this is a strong function of planet mass and 
radius (as the transition occurs when R p /R s -C 1). Sim¬ 
ilarly, close-in low mass planets will begin evaporating in 
the recombination limit as they contract significantly due 
to cooling and the flux drops they may make a transition 
to photon-limited evaporation; however as the flux drops 
further and they contract even more they will eventually 
transition to energy limited evaporation. Whether or not 
low-mass planets experience photon-limited evaporation 
and when they make the transition from recombination 
to photon to energy limited is non-trivial. What is cer¬ 
tain is there is not a unique value of the flux that one 
transitions from one regime to the other as assumed in 
previous evolutionary studies (e.g. Jin et al. 2014). Fur¬ 
thermore, our analysis presented here makes no effort to 
quantify what the efficiency of energy limited flow is and 
that can only be answered by numerical calculations such 
as those we preform in the next Section. 

3. NUMERICAL SIMULATIONS 

In the previous section we have put forward the basic 
physical picture of planetary evaporation driven by EUV 
radiation, and their scalings. However, the actual mass- 
loss rates are still unknown as the details of the picture 
remain unquantified. In particular, the efficiency of any 
energy limited flow, or the exact temperature profile of 
any radiation recombination limited flow is sensitive to 
the exact details of the heating and cooling. 

Furthermore, since we only expect the ionization front 
to be thin in the case of a flow which is recombination 
limited, and as we shall see it often the scale of the planet, 
one cannot assume a monochromatic irradiation field (as 
done by Murray-Clay et al. 2009) and radiation harden¬ 
ing will need to be taken into account. Thus, the only 
way to obtain an accurate picture of planetary evapo¬ 
ration we must proceed and perform numerical calcula¬ 
tions. 











6 


Owen, J. E. & Alvarez, M. A. 






Fig. 1.— The various mass-loss regions are shown as a function of planet radius and mass for various incident EUV fluxes and mass-loss 
efficiencies. The crossed green region shows “photon-limited” mass-loss, the blue striped region show “energy-limited” mass-loss and the 
remaing parameter space is recombination limited. The thin lines show lines of constant density at 0.1, 1.0 & 10 g cm _l! and the points 
show currently observed exoplanets obtained from the open exoplanet catalogue - data obtained on 5th February 2015 - (Rein, 2012). The 
top left panel shows our fiducial case for a flux of 3 X 10 4 erg s -1 cm -2 and efficiency r/ = 0.1. The bottom left and top right panels show 
fluxes of 3 X 10 5 & 3 X 10 3 erg s —1 cm -2 respectively. The bottom right panel shows a flux of 3 X 10 4 erg s —1 cm -2 and efficiency // = 0.3. 


We follow many other previous approaches and con¬ 
sider ID spherically symmetric outflow, where we solve 
the problem along a streamline connecting the star 
and planet, neglecting the Coriolis force. Along this 
streamline we solve the radiation-hydrodynamic problem 
including EUV irradiation, and non-equilibrium heat¬ 
ing/cooling and Hydrogen ionization chemistry. We use 
the ZEUS-MP hydrodynamics code (Stone & Norman, 
1992, Hayes et al. 2006), which has been previously used 
for planetary evaporation studies (Owen & Adams, 2014) 
and add to it a new module to perform the radiative- 
transfer, ionization structure and heating and cooling 

update. _ 

As discussed in Section |2.1.1| our numerical setup will 
result in a “photon-limited” mass-loss rate that is too 
high compared to reality. In fact in all simulations we 
perform we find that the photon limited rate of our spher¬ 
ically symmetric setup is much higher than any measured 
mass-loss rates. As such in our numerical study here we 
only investigate the transition between energy and re¬ 
combination limited and we must leave an investigation 


of the transition to photon limited flows to future multi¬ 
dimensional calculations. 

3.1. Ionizing radiative transfer 

Ionizing radiative transfer can be solved most effi¬ 
ciently in a causal manner where one performs the ray 
tracing calculation proceeding away from the source (to¬ 
wards the planet in our case) in a step-wise manner. We 
operator split our radiative transfer calculation from the 
hydrodynamics and perform it between the “source” and 
“transport” steps of the standard ZEUS scheme. We con¬ 
sider only pure Hydrogen and the equations for ionization 
and energy due to this update are given by: 

-—= (1 — X)(Fi + n e Cn) — Xn e aB ( 21 ) 

at 

dc 

= -PX7 . U + T h -A (22) 
at 

where I\ is the photoionization rate, Ch is the colli- 
sional ionization rate, e is the internal energy of the gas, 












UV driven evaporation of planets 


7 


T/j is the heating rate and A is the cooling rate. We 
include heating from photo-ionizations and cooling from 
recombination, collisional excitation, collisional ioniza¬ 
tion and free-free emission (Katz et al. 1996). How¬ 
ever, we neglect heat transport due to conduction as it 
known to be sub-dominant by many orders of magnitude 
(e.g. Murray-Clay et al. 2009). We follow (Mellema 
et al. 2006) and pre-tabulate the heating and ioniza¬ 
tion rates as a function of neutral hydrogen column den¬ 
sity for our input spectrum, thereby capturing the non- 
monochromatic nature of our problem. Equations [21] & 
[22] are then solved explicitly by sub-cycling then on time- 
scales which can be shorter than the hydrodynamic time- 
step. We note that our scheme requires that we resolve 
the ionization front, something we have enough resolu¬ 
tion to do. 

3.2. Numerical setup 

We perform a series of simulations in order to bet¬ 
ter understand the regimes which planetary evaporation 
exists in. Since a full exploration of the entire parame¬ 
ter space numerically is beyond the scope of this work, 
we instead choose to focus on low-mass planets (where 
evaporation is more evolutionary important - see Owen 
& Wu 2013; Lopez & Fortney 2013). In all simulations 
we use an underlying bolometrically heated atmosphere 
of 1000 K, which is forced to be isothermal by means 
of a temperature floor. We adopt a grid with 192 non- 
uniformly spaced cells with a resolution of ~20 km at 
the base of the grid, and a radial extent from R p to 
30 R p . The simulation is initiated with a neutral hy¬ 
drostatic density structure for a 1000 K atmosphere. At 
the lower boundary the ghost cells are reset to the neu¬ 
tral, hydrostatic density and temperature structure ev¬ 
ery time-step. At the upper boundary we use standard 
outflow boundary conditions, which are perfect in ID 
super-sonic outflow (see Stone & Norman 1992; Hayes 
et al. 2006, for their implementation in the ZEUS code). 
In order to choose the density at the inner boundary 
we set the density, such that the initial neutral hydro¬ 
static atmosphere has an optical depth of r = 100 to 
13.6eV photons, which typically results in densities of or¬ 
der ~ 10 -11 — 10~ 12 g cm -3 . We find that such a choice is 
sufficient to keep the inner regions of the active domain 
part of the bolometrically heated atmosphere rather than 
UV heated atmosphere. The domain is then irradiated 
by the UV and a transonic thermal wind is launched. In 
all cases the flow approaches steady-state and no vari¬ 
ability is found. All results are analysed after ~ 100 
sound crossing times after steady-state is firmly estab¬ 
lished, with the mass-flux constant at the ~ 10~ 3 level 
(see Figure 2). 

We perform two series of simulations: one at con¬ 
stant planet mass with a mass of 10 M® where we 
vary the planet radius and one set at constant radius of 
3.14 R® where we vary the planet’s mass. In both cases 
the radius/mass are varied such that the smallest/most 
massive planets have a density of 10 g cm~ 3 and the 
largest/least massive have at ratio of planet radius to 
the “Parker” radius^] of 0.2. For each set we perform 24 
for individual simulations at 12 different fluxes varying 

3 The “Parker” radius is the sonic radius of an isothermal flow, 
where we take the temperature of the underlying bolometrically 


between ~ 10 and ~ 10 6 erg s” 1 cm -2 (where we take 
the EUV luminosity contained between 13.6 and 30eV). 
Our input spectrum was calculated by Ercolano et al. 
(2009) and has been used in previous exoplanet evapora¬ 
tion studies by Owen & Jackson (2012) and extends out 
to several lOOeV (although photons with high energies 
have an incredibly small photo-ionization cross-section 
with Hydrogen, a oc E~ 3 ). 

3.3. Results 

As expected we identify two limiting flow types: re¬ 
combination limited at high fluxes and energy limited 
at low fluxes. Note as discussed above “photon-limited” 
flow cannot be found for realistic parameters in ID spher¬ 
ically symmetric simulations. We show the different flow 
topologies for recombination limited (top panels) and en¬ 
ergy limited flow (bottom panels) in Figure [2] Both sim¬ 
ulations are for a planet with a radius of 3T4 R® and 
mass of 6.25 M®. The recombination limited flow has a 
flux of 10 6 erg s _1 cm -2 and the energy limited flow has 
a flux of 10 3 erg s _1 cm -2 . 

In the “recombination-limited” case (top panels), we 
see that the size of the ionization front is small compared 
to the planet’s radius and the flow is fully ionized at large 
radius indicating recombination-ionization equilibrium. 
The temperature quickly reaches 10 4 K at the ionization 
front and the sonic point is very close to the planet’s 
surface. Furthermore, we note that most of the heating 
is concentrated close to the planet’s surface, as indicated 
by the heating rate profile ( qr 3 - multiplied by r 3 to 
indicate the radius where the majority of heating takes 
place) shown in the bottom right panel of Figure 2a. We 
find that at these low masses the flow does adiabatically 
cool slightly at large radius although this is in the super¬ 
sonic regime of the flow it doesn’t affect the mass-loss 
rate when compared to the simple estimates presented 
above. 

In the energy limited case (bottom panels), we find 
that the ionization front essentially encompasses the en¬ 
tire flow and the flow is still partly neutral at 10 planetary 
radii. This results in a much colder flow only reaching 
~ 5000 K at 2 planetary radii, before adiabatically cool¬ 
ing at large radii. This results in the sonic point being 
pushed to a larger radius. This flow topology is empha¬ 
sised by the fact the heating rate peaks at larger radii 
than in the recombination limited case and is more spa¬ 
tially extended, as shown in the bottom right panel of 
Figure 2b. 

We can investigate the transition between energy- 
limited and recombination limited in Figure [3] where we 
plot the normalised mass-loss rate (where the mass-loss 
rate for all models is normalised such the mass-loss rate 
is unity at the lowest flux of ~ 10 erg s _1 cm -2 ) against 
flux (top panel) and mass-loss rate divided by flux (bot¬ 
tom panel) for a planet with several different masses: 
5.12,6.25, 10.31, 13.92 M® (with a radius of 3.14 R®). In 
the insert we show the logarithmic derivative of mass-loss 
rate with flux. At low fluxes the mass-loss rate scales lin¬ 
early with flux following the dashed line as expected for 
energy limited flow; however, at large fluxes the mass-loss 

heated atmosphere and is given by: GM P / 2c^; planets larger than 
this radius are not hydrostatically bound over evolutionary time- 
scales. 
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(b) 

Fig. 2.— Flow toplogy for a planet with radius of 3.14 R ,0 and mass of 6.25 M 0 receiving an ionizing flux of 10 6 erg s -1 cm -2 (a) and 
10 3 erg s -1 cm -2 (b). The top left panels show the gas temperature; bottom left panels show the ionization fraction; top middle panels 
show the velocity (solid) and sound speed (dashed); bottom middle panels show the mass-flux (i.e. 47 rr 2 pu); top right panels show the 
density and bottom right panels show the volumetric heating rate ( q ) multiplied by r 3 to indicate where most of the heating is taking place. 
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Fig. 3.— The top panel shows the normalised mass-loss rate 
(the mass-loss rate for all cases is scaled to equal unity at the 
lowest flux of ~10 erg s -1 cm -2 ) as a function of flux shown as 
solid lines for a planet with a radius of 3.14 Rq and masses of 
(from top to bottom) 5.12, 6.25, 10.31, 13.92 M®. The dashed line 
shows the scaling expected for a pure “energy” limited flow. The 
dotted line show lines propotional to the square root of flux, the 
scaling expected from recombination limited flows. The bottom 
panel shows the same mass-loss rates divided by flux (again scaled 
so that the ratio of mass-loss rate and flux is one at the lowest flux 
of ~10 erg s -1 cm -2 . The insert shows the logarthmic derivative 
of mass-loss rate by flux at the different masses 5.12 (solid), 6.25 
(dashed), 10.31 (dotted), 13.92 (dashed) M®, the vertical lines 
show the analytical parscription would predict the flow transitions 
from energy limited to recombination limited. 

rate scales approximately as the square root of flux (fol¬ 
lowing the dotted lines), indicating recombination lim¬ 
ited flow. We find that the transition point varies with 
mass and perhaps more interestingly the transition is 
slow, taking a large range in flux. The predicted tran¬ 
sition fluxes from Section 2 - as indicated by the vertical 
lines in the insert of Figure [3] - pick out the end of the 
transition. As expected planets with higher masses be¬ 
come recombination limited at lower fluxes due to there 
increased flow-time scales. 

We can see this transition in more detail if we plot 
the mass-loss rate as a function of flux for various planet 
masses, for the full simulated range of incident flux. The 
results are shown in Figure [4] with the 1/M p scaling 
for energy-limited flows shown as dotted lines. We find 
at low fluxes the simulation is well matched by an en¬ 
ergy limited scaling; however, at higher fluxes the trend 
becomes much flatter indicating recombination limited 
flows (Noting that in this parameter range the planet’s 
escape temperature is < 10 4 K so mass plays a limited 
role as discussed above). The dashed line in Figure [4] 
shows the predicted analytic transition from recombina¬ 
tion to energy limited flows. Above the line is predicted 
to be recombination limited and below is predicted to be 
energy limited. We see this prediction agrees well with 



Fig. 4. — Mass-loss rate as a function of planet mass for the full 
ranges of fluxes (bottom to top) logarithmically spaced between 10 
and 10 6 erg s -1 cm -2 . The dotted lines show the 1 /M p scaling 
expected from “energy” limited flows. The dashed line shows the 
predicted transition between recombination limited flows and en¬ 
ergy limited flows from Section 2, above the line is predcited to be 
recombination limited and below is predicted to be energy limited. 



Fig. 5. — Mass-loss rate as a function of planet radius. Simulation 
results are shown as points for fluxes (bottom to top) of ~ 10, 30, 
100 erg s —1 cm -2 . The solid lines shows the “energy” limited 
mass-loss rate with an efficiency of 20%. 

the numerical results. 

Finally, in order to measure the efficiency of the flow 
we measure the mass-loss rate as a function of planet 
radius (due to its steep dependence this should provide 
a good range of mass-loss rates over which to tune this 
parameter). The discussion in Section 2 indicates that 
at low fluxes and low masses the flow should be strongly 
energy limited, which is indeed what we find. In Figure [5] 
we show the mass-loss rate as a function of radius for a 
planet of mass 10 M® and fluxes of ~ 10-100 (points). 
We find that an mass-loss rate with an energy efficiency 
of 20% agrees well with the simulation^ and as such 

4 Note in order to compare with Equation [5] we must drop the 
factor 4 in the denominator that arises from geometrical averaging 
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represents a suitable measure of the mass-loss efficiency 
in the region of parameter space. 

4. DISCUSSION 

We have argued that evaporation due to UV radiation 
can exist in three limiting cases: “recombination-limited” 
- where the flow is in radiative-recombination equilibrium 
and approximately isothermal at 10 4 K; “energy-limited” 
where incoming radiation is balanced by mechanical lu¬ 
minosity; and a new regime: “photon-limited” - where the 
mass-loss is limited by the fact that each ionizing photon 
can allow one hydrogen atom to escape. Each distinctive 
limit has a different scaling with planet mass, radius as 
well as the incident flux. We have shown that there is no 
unique transition flux from one regime to the other and 
such a limit depends on the planetary parameters. 

Schematically if the flux is high enough such the the 
density in the flow is large (and correspondingly the re¬ 
combination time is short) or the planet is so big or mas¬ 
sive that the flow time-scale is long then the problem is 
in the radiative-recombination limited, where the recom¬ 
bination time is much faster than the flow-time-scale. As 
argued by Murray-Clay et al. (2009); Owen & Jackson 
(2012) in such a scenario it behaves as a HII region that 
is approximately isothermal and the density at the ion¬ 
ization front scales as J ^ 2 , in this region the mass-loss 
rate (in a ID calculation) is analytically tractable and 
we present such a calculation in Section 2.1.2. 

If one decreases the flux, the density in the ionization 
front will drop and the recombination time will increase. 
Once it reaches the point where the recombination time 
is longer than the flow time-scale then recombinations 
can be ignored. When the recombination time is much 
longer than the flow time then we also find the cooling 
time is longer than the recombination time and subse¬ 
quently the flow time-scale. In this case radiative cooling 
can be ignored. This is what one thinks off as standard 
“energy-limited” flow. However, we note that when the 
recombination time is of order the flow-time scale then 
the cooling time is shorter than the recombination time 
and thus shorter than the flow time-scale. Thus, one 
expects the transition from recombination to energy lim¬ 
ited to happen smoothly and broadly over a large range 
of fluxes. This is indeed what our numerical calculations 
indicate, that the transition occurs over several orders of 
magnitude in flux. This is different to the case of X-ray 
heating discussed by Owen & Jackson (2012), where since 
the heating and cooling times scaled identically with flux 
then the cooling time was always comparable or shorter 
than the flow time-scale. Thus, in the case of X-ray heat¬ 
ing Owen & Jackson (2012) argued that the flow was 
never “energy-limited”. 

If when one reaches the point where the recombination 
time is becoming longer than the flow time-scale and the 
planet is not tightly bound (such that the energy limit 
is not restrictive) then the mass-loss rate will be limited 
by the number of ionizing photons. Such that one can 
only lose one hydrogen atom per ionizing photon. Obvi¬ 
ously, in this limit planet mass is completely irrelevant 
and the mass-loss rate just scales as the flux x area. 
Such a case is similar to the evaporation of low-mass 
neutral clumps by nearby young massive stars discussed 

which does not occur in spherically symmetric simulations. 


by Bertoldi & Mckee (1990). However, this regime will 
only be restricted to a small (but perhaps important) 
regime of parameter space for low-mass large radii plan¬ 
ets, similar to the masses and radii of the Kepler planets 
at young 10s Myr ages (e.g. Lopez et al. 2012, Owen 
& Wu, 2013). Unfortunately, such a problem could not 
be investigated with realistic planet parameters in ID 
spherically symmetric simulations. This is because such 
simulations assume spherical irradiation, in this case the 
absorbing area would be 47ri?| F . However, in the reality 
of planet-parallel irradiation and the fact the ionization 
front is marginally optically thin then we suspect the ab¬ 
sorbing area is more likely ~ nR^. In this regime the 
ionization front is thick and by construction has a radius 
i?iF R p and as such the “photon-limited” mass-loss 
rate from the simulations is much higher than reality and 
any other limit. One needs to perform 2/3D simulations 
to investigate the role of “photon-limited” mass-loss in 
planetary mass-loss. 

Our numerical calculations indicate that the hydro- 
dynamic efficiency if fairly low for low-mass planets 
~ 10 — 20 %, comparable to the values used in plane¬ 
tary evolution studies (e.g. Lopez et al. 2012, Lopez & 
Fortney 2013) indicating that much of the absorbed en¬ 
ergy remains as kinetic/thermal energy in the flow rather 
than being converted to mass-loss. When the flow tran¬ 
sitions to energy or photon limited the ionization front 
is thick; essentially encompassing the entire the flow. 
Furthermore, we note that not all the heating is con¬ 
centrated at the planet’s surface as can be seen in the 
bottom panels of Figure [2j where the maximum tem¬ 
perature is reached at ~ 2 planetary radii. The ioniza¬ 
tion fraction increases monotonically with radius but still 
contains a significant neutral fraction out to many plan¬ 
etary radii. Therefore simplified calculations assuming 
infinitesimally thin ionization fronts are not appropriate 
in this limit. This means photon hardening can play an 
important role in pre-heating material; this occurs when 
harder photons penetrate deeper than the standard ion¬ 
izing (~ 13.6 eV) photons resulting in an average photon 
energy that is harder than the total spectrum, these pho¬ 
tons ’pre-heat’ the underlying gas. Thus, one requires a 
non-monchromatic numerical method to take this into 
account (unlike Murray-Clay et al. 2009, who assume a 
spectrum with just 20 eV photons). 

Finally, we note we have presented a new numerical 
algorithm for the inclusion of UV heating in planetary 
evaporation simulations. We have included the frequency 
dependence of the ionizing spectrum for pure Hydrogen 
atmospheres and note this can easily be extended to in¬ 
clude more species (Friedrich et al. 2012). Such an al¬ 
gorithm can be incorporated into multi-dimensional cal¬ 
culations that include ray-tracing (e.g. Owen & Adams, 
2014) to investigate planetary evaporation in 2/3D. Such 
a further step is important, not only to investigate the 
flow properties, in particular “photon-limited” evapora¬ 
tion identified here; but, also to produce realistic calcu¬ 
lations that can be compared to observational indicators 
of planetary evaporation. This includes measured tran¬ 
sits in Lyman-a (Vidal-Madjar et al. 2003, Ben-Jaffel et 
al. 2007, 2008, Lecavelier Des Etangs et al. 2010) and 
other metal lines (Vidal-Madjar et al. 2004, Ballester & 
Ben-Jaffel 2015). 
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5. SUMMARY 

We have investigate planetary evaporation under the 
influence of irradiation by ionizing UV photons at close 
enough separations such that evaporation takes place in 
the hydrodynamic limit. We argue that there are three 
main evaporation limits. The first two were demon¬ 
strated by Murray-Clay et al. (2009), at high fluxes the 
flow is in radiative-recombination equilibrium and the 
mass-loss rate scales as the square-root of the ionizing 
flux; the second where recombinations are negligible and 
the mass-loss rate is limited by the amount of incom¬ 
ing energy and as such scales linearly with the ionizing 
flux. Finally, we argue that there is a third limiting factor 
that must be taken into account, which we term “photon- 
limited” mass-loss. In photon-limited mass-loss one can 
only lose a single hydrogen atom per incoming ionizing 
photon and is important when the planets gravity is weak 
(low density, low mass planets). 

We present a basic physical picture of the three differ¬ 
ent evaporation regimes and show that transitions from 
one regime to another does not happen a specific UV 
flux, as assumed in some recent evolutionary studies (e.g. 
Jin et al. 2014). The different evaporation regimes domi¬ 
nate different planetary parameters and different flux and 
we present a framework that can be incorporated in to 
planet evolution studies, provided a mass-loss efficiency 
is chosen for the energy-limited flows. 

We have presented a new numerical algorithm for per¬ 
forming radiation-hydrodynamic simulations that takes 
into account the frequency dependence of the irradiation 
spectrum. It does not require the ionization front to be 
thin, or approximately isothermal and self consistently 
finds the correct ionization and thermodynamic struc¬ 
ture of the flow. This radiative transfer scheme has been 
incorporated into the ZEUS hydrodynamic code. 

Using ID spherically symmetric numerical simulations 
we have investigated the transition between “energy- 
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limited” flows and “recombination-limited” flows. We 
take into account the non-monochromatic nature of the 
ionizing photons unlike previous works (e.g. Murray- 
Clay et al. 2009). The transition from recombination 
limited to energy limited occurs when the recombination 
time becomes longer than the flow time-scale and the 
ionization front encompasses the entire flow scale. In 
this case radiation hardening from the full spectrum is 
an important contribution to the thermodynamics. 

We find that the transition from recombination lim¬ 
ited to energy limited is broad and takes several orders of 
magnitude in flux. Planets with deeper potentials tran¬ 
sition to energy-limited flows at higher fluxes consistent 
with our analytically expectations. Furthermore, we find 
that in the‘energy-limited case efficiencies around 10-20% 
provide a good fit to mass-loss rates. 

Finally, we could not use our ID spherically symmet¬ 
ric simulations to investigate “photon-limited” flows and 
future mutli-dimensional calculations are required to in¬ 
vestigate the importance of this type of flow. Something 
that is possible with our presented algorithm and we plan 
to incorporate this into the multi-dimensional ray-tracing 
scheme for planetary evaporation presented by Owen & 
Adams, (2014). 
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